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We explore the interplay of membrane curvature and non- 
specific binding due to excluded-volume effects among col- 
loidal particles inside lipid bilayer vesicles. We trapped sub- 
micron spheres of two different sizes inside a pear-shaped, 
multilamellar vesicle and found the larger spheres to be 
pinned to the vesicle's surface and pushed in the direction 
of increasing curvature. A simple model predicts that hard 
spheres can induce shape changes in flexible vesicles. The 
results demonstrate an important relationship between the 
shape of a vesicle or pore and the arrangement of particles 
within it. 

Entropic cxcludcd-volumc (depletion) effects are well 
known to lead to phase separation in the bulk of colloids 
and emulsions consisting of large and small particles with 
short-range repulsive interactions More recently, 

attraction of the large particles to flat, hard walls [QJ|] 
and repulsion from step edges || have been demonstrated 
in binary hard-sphere mixtures. A key concept suggested 
in these papers is that the geometric features of the sur- 
face can create "entropic force fields" that trap, repel or 
induce drift of the larger particles. This mechanism is 
not limited to suspensions of micron-sized particles; it 
may play a role in "lock and key" steric interactions on 
smaller macromolecular length scales. For example, the 
shape of pores and liposomes inside cells is likely to af- 
fect the behavior of macromolecules confined within them 




In this Letter, we present experimental results that 
demonstrate new entropic effects at surfaces. In particu- 
lar, the behavior of particles confined within vesicles re- 
veals quantitatively the striking effect of membrane cur- 
vature. We first discuss experiments probing the behav- 
ior of a microscopic sphere trapped inside a rigid, phos- 
pholipid vesicle. Adding much smaller spheres to the 
mixture changes the distribution of the larger sphere in a 
way that depends on the curvature of the vesicle wall (see 
Fig. |](b) and (c)). The results are consistent with the 
depletion force theory and illustrate a new mechanism for 
the size-dependent arrangement of particles within pores. 
We then explore theoretically some consequences of re- 
placing the rigid wall with a flexible one. The entropic 
curvature effects can overcome the membrane's stiffness, 
leading to a new mechanism for shape changes in vesicles. 

We first briefly review depletion effects in mixtures of 



microscopic hard spheres of two different sizes. Moving 
two of the larger spheres toward one another does not 
change their interaction energy (which is zero for hard 
spheres) but does increase the volume accessible to the 
other particles (Fig. ^|). The resulting gain in entropy 
reduces the free energy of the system by (^/2)a<j)sk^T 
flUp^l . Here, a is the ratio of large to small radii 
(Rl/Rs), 05 is the small-sphere volume fraction, and 
ksT is Boltzmann's constant times the absolute temper- 
ature. This simple result relies on the approximation that 
the small spheres are a structureless ideal gas and that 
the large-sphere volume fraction, </>£, is small. The re- 
duction of free energy produces an "entropic force" that 
pushes the large spheres together. When the large sphere 
is moved to a flat wall, moreover, the overlap volume and 
the free-energy loss are approximately doubled . In bi- 
nary hard-sphere mixtures, these effects are known to 
drive crystallization of large spheres in the bulk jl],|]||] 
and at a flat surface |I],|7],|l3],|l4| ■ Furthermore, the shape 
of the wall can lead to entropic forces in a specific direc- 
tion along the wall. For example, the larger spheres are 
locally repelled from an edge cut into the wall [gj and 
attracted to a corner (i.e. where the "wall" meets the 
"floor") |l5| |. If the wall has constantly-changing radius 
of curvature, these forces are predicted to act everywhere 
along it Q. As shown in Fig. || (c), when the large sphere 
is near the wall, the overlap volume depends on the wall's 
curvature radius. The large sphere will therefore move 
in the direction of increasing curvature to minimize the 
small spheres' excluded volume. 

To measure this surface entropic force, we have studied 
the distribution and dynamics of microscopic nearly-hard 
spheres trapped inside rigid, pear-shaped vesicles. The 
vesicles were prepared from a phospholipid (1-Stearoyl- 
2-01eoyl-sn-Glycero-3-Phosphocholine (SOPC), Avanti 
Inc., U.S.A.), dissolved in chloroform (25 mg/mL). Af- 
ter evaporating the chloroform from 200 /iL of SOPC 
solution, we added 100 [iL of salt water with charge- 
stabilized polystyrene spheres (Seradyn, IN, USA) in sus- 
pension. The salt (0.01 M NaCl) served to screen out 
electrostatic forces over a distance of « 5 nm, the Debye- 
Huckel screening length [Q. Thus, our sample closely re- 
sembled an ideal hard-sphere and hard-wall (HSHW) sys- 
tem. Rigid, multilamellar vesicles of diverse shapes and 
sizes immediately formed with colloidal spheres trapped 
inside. We injected the solution into a 10-/xm thick glass 
container for viewing under an optical microscope (lOOx 
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objective with 1.30 numerical aperture, in transmission 
mode). Images of a planar slice, 600 nm thick, through 
the center of vesicles were captured, then later digitized 
(Fig. 0(a)). _ 

We quantified the behavior of the mixture by mea- 
suring its free energy, F(r), as a function of the posi- 
tion of a sphere of radius Rl = 0.237/im. First, using 
NIH Image software, we determined the in-plane posi- 
tion of the sphere's center of mass, r, (±0.08/mi) when 
it appeared in the imaged slice. From these data, we 
used two techniques to extract F(r). In the first, the 
number of times Nfa) the sphere appeared at a given 
bin located at rj defined the density distribution. As- 
suming that each measurement event was independent 
of the others, N(ri) follows the Boltzmann distribution: 
iV(rj) oc exp(— F^r^/k^T). Since a systematic error can 
arise if events are not completely independent (i.e. not 
separated by an infinite time), we waited a minimum of 
0.6 s between measurements. During this time the mean 
square displacement was 0.3 /im, larger than the 0.07- 
0.13 /mi bin sizes. We also collected data over a period of 
30-80 minutes, enough time for the sphere to explore all 
of the available space. We therefore used the logarithm 
of the sphere's distribution to extract the free energy. 

The dynamics of the diffusing 0.237-/im sphere pro- 
vided the second way to measure F(ri), as described in 
Jl6| . The region along the inner surface of the vesicle was 
divided into several equal-area bins. We considered only 
events in which the center of the large sphere was within 
0.28 /im of contact with the surface in consecutive frames 
(separated by a time r). From our videotape, we counted 
the number of times (Nij(r)) the sphere was located in 
bin j at time t and in bin i at time t + r. The transi- 
tion probability matrix, Pjj(r), is given by Nij{r)/N(rj). 
The measured matrix contains information about the 
equilibrium state of the system. In particular, the eigen- 
vector (e) of Pij with unit eigenvalue is proportional 
to the Boltzmann factor: e\; oc exp(— F(ri)/k-QT). We 
obtained consistent results with r = 0.1, 0.2 and 0.3 s 
and with various bin sizes. This technique avoids fit pa- 
rameters and possible systematic errors arising from the 
density-distribution approach. It also avoids potential 
errors arising from a slow change in shape of the vesi- 
cle, which would affect the distribution averaged over a 
long time but would not affect the short-time particle 
dynamics. We report here the results from two different 
samples. The first "control" sample contained a solitary 
0.237-/im sphere (no small spheres) diffusing freely inside 
a vesicle. The measured density distribution is shown in 
Fig. 0(b). The sphere distribution is uniform: there is no 
significant interaction between the vesicle wall and the 
polystyrene sphere, as expected due to the very short 
Debye-Huckel screening length. 

The behavior changed noticeably when small spheres 
were added to the interior of a vesicle. In Fig. 0(c), 
we show the distribution of the Rl = 0.237 /im sphere 



in a binary mixture with (f>s = 0.3, Rs — 0.042/iin 
(a = 5.7). The distribution is highly non-uniform, with a 
significantly higher probability of finding the large sphere 
within about 0.28 /tm of the surface. The apparent width 
of this "surface" region exceeds 2R$ due to uncertainties 
in the particle positions and, possibly, due to a slight tilt 
of the vesicle wall away from vertical. 

We measured the average number of times the large 
sphere appeared in each bin within 0.28 /tm of the sur- 
face and the average number per bin in the bulk. We 
defined the natural logarithm of the ratio of these num- 
bers as F /(kBT), a measure of the average strength of 
the depletion attraction. We found F = (2.2 ± 0.5)fcsT. 

Theoretically, 

F = In J^ +0 ' 28 Mm [dr/0.28 fim] e^MAsT where r is 
the distance from the large-sphere center to the wall and 
V(r) is given by the depletion force model. Although the 
vesicle wall was curved, we can predict a lower bound of 
Fq by assuming the wall is flat: 
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R L -2R s ) 2 (r + 2R L + R s ). (1) 



Putting in the numbers gives a theoretical prediction of 
Fq = 1. 97 ksT. Here we have neglected thermal fluctu- 
ations in the shape of the vesicle wall and residual elec- 
trostatic repulsion. Constraining the large sphere to lie 
5 nm away from the wall, for example, reduces the pre- 
dicted result to 1.45 k B T. The 0.7-fc B T difference from 
our measured result is likely due to the curvature of the 
wall, which would enhance the observed Fq as discussed 
in the following paragraph. These results are consistent 
with recent calculations showing that the ideal-gas ap- 
proximation accurately predicts the depletion well depth 
at contact (although it misses the relatively weak, long- 
range depletion repulsion) |p_ 7|| . 

The distribution in Fig. [j[(c) demonstrates a higher 
large-sphere occupation where the vesicle surface is more 
curved. ^From the sphere dynamics, we obtained the 
total free energy, F, as a function of sphere position 
when it was near the surface. The inset of Fig. || shows 
the measured F as a function of the position, s, along 
the perimeter. The slope of this curve reveals a maxi- 
mum force of 20xlO~ 15 iV pushing the large sphere along 
the wall. We also determined the curvature radius at 
points along the wall using our images of the vesicle. 
In Fig. we plot the free energy as a function of wall 
curvature radius, Rc- The free energy decreased by 
approximately 1.5fceT when the curvature radius de- 
creased from 20 /tm to 2 /im. For comparison to our 
results, we used the simple depletion-force model de- 
scribed above (Fig. ||(c)) to calculate the free energy as a 
function of wall curvature (assuming a locally spherical 
wall): F(R C ) = -Ps Overlap (-Rc). Here, P s is the small- 
particle osmotic pressure (from the ideal-gas law) and 
Overlap {Rc) is the size of the excluded- volume overlap. 
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The large sphere was assumed to touch the vesicle wall. 
The theory, represented by the solid line in Fig. ^, is con- 
sistent with the experimental results. For both plots, the 
free energy was set to in the limit of large curvature 
radius (flat wall). Since we could only measure Rc in 
the image plane, we assumed that this also determined 
the out-of-plane curvature. This approximation probably 
explains the scatter in our data points. 

So far we have shown that curvature of a rigid wall in- 
duces forces on the large spheres in binary particle mix- 
tures. To conclude this Letter, we consider theoretically 
a binary hard-sphere mixture in the presence of flexible 
walls. This extra degree of freedom introduces a compe- 
tition between depletion and curvature energy that can 
produce a variety of new phenomena. For example, we 
expect that under some circumstances, the membrane 
will spontaneously bend around the large sphere (Fig. 
We anticipate that these effects will be observable and 
that understanding this mechanism will be essential to 
predict the behavior of unilamellar phospholipid vesicles 
containing particles of different sizes. 

The membrane will spontaneously envelop the large 
sphere only if the resulting total free energy change, 
AF, is negative. Our experiment indicates that the 
membrane-particle interactions are very short-ranged, so 
we use the hard-sphere and hard-wall approximations. 
We therefore divide the full free energy into two terms, 
one representing the net adhesion induced by the col- 
loidal spheres, the other the vesicle's bending energy: 

R -^adh "1" Fves • 

The experiment described above shows that for rigid 
membranes, -F a dh can be well described by the depletion- 
force model. For flexible membranes we must add a new 
contribution to -Fadh to account for the "steric repulsion" 
effect For R$ -C Rl we may approximate the lat- 

ter by the interaction between a flat hard wall and a fluc- 
tuating membrane at constant mean separation x. This 
free energy equals 0.06(/cbT 1 ) 2 / ' nx 2 per unit area po| , ^l| . 
For unilamellar phospholipid membranes, the bending 
stiffness k ~ 15/sbT |22]. Approximating the change of 
excluded volume upon adhesion by 2R$ — x per unit area 
and again taking the small-sphere volume fraction to be 
4>s = 0.3, we minimize the total adhesive energy (deple- 
tion attraction plus steric repulsion) to find x ~ 0.5i?s 
and the adhesive energy F adh = -0.09lk B T ■ A*/R s 2 , 
where A* is the area of contact between the membrane 
and the sphere (Fig. |(b)) ||. 

The vesicle's elastic energy, F vcs , depends only on the 
membrane curvature. We neglect the constant-volume 
and constant-area constraints of a vesicle and instead 
consider an infinite membrane. Thus the local curvature 
energy change upon adhesion is just F vcs = 2kA*/Rl 2 
p3,p5l. Combining F a( jh + -Fvcs, we see that adhe- 
sion and engulfment are favored when a = Rl/Rs > 
■v/2K/0.091fcBT =18. A more detailed analysis of the 
deformation of vesicles at a generically sticky surface, 



specialized to the case of interest here, gives a similar re- 
sult [ p6[ . The required regime for a is easily accessible, 
using either polymers or spheres as the small particles. 

In conclusion, we have measured for the first time 
a curvature-induced entropic force in a system of hard 
spheres trapped in a rigid vesicle. The results show that 
the distribution of particles within the vesicle is strongly 
affected by the local shape of the vesicle wall. Further- 
more, a simple estimate predicts shape transitions of unil- 
amellar phospholipid vesicles induced by particles inside 
the vesicle. The ideas presented here suggest a way to 
understand several phenomena in cellular interiors and 
complex fluids inside porous media. For example, in 
a vesicle whose membrane contains multiple species of 
lipids, the lipids can segregate into regions with different 
curvatures p7| . In such a sample, the particle distribu- 
tion, which depends on curvature, would correlate with 
the local composition of the membrane. This mechanism 
may provide considerable specificity in associations be- 
tween macromolecules and lipids. Finally, when the size 
ratio is more extreme or the vesicle wall is less rigid, the 
spheres may be able to induce budding and fission of the 
vesicle. 
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FIG. 1. (a) Cartoon of the 600-nm thick slice through an SOPC vesicle imaged with an optical microscope. We measured 
the in-plane positions of the larger colloidal sphere when it was in focus, (b) (color) Probability distribution of a single 
0.237-^tm-radius polystyrene sphere inside a vesicle (no small spheres). The white dashed line is the edge of the vesicle, and 
the colored points indicate the number of times, N, the center of the sphere was observed in a bin located at a given point. 
There were 2000 events and the bins were 130 nm x 130 nm. The sphere simply diffused freely throughout all of the available 
space, (c) (color) Same as in (b), but with a vesicle that also contained small spheres (4>s = 0.30, Rs = 0.042 /im,a = 5.7). 
There were 2300 events and the bin size was 65 nm. The large sphere was clearly attracted to the vesicle wall, especially where 
the vesicle was most curved. 
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FIG. 2. Description of depletion forces in a binary 
hard-sphere mixture. The centers of mass of the small spheres 
are excluded from the hatched regions, within one radius of 
the surface of the large sphere and of the walls. In (a), the 
volume accessible to the small-sphere centers, V a cc is the total 
volume minus the hatched regions, (b) When the large sphere 
moves to the wall, the excluded regions overlap, as shown in 
black, and 14cc increases by this amount. The small spheres' 
entropy therefore increases, (c) Because of the changing wall 
curvature, the size of the overlap depends on sphere position. 
The large sphere will move along the wall to maximize the 
size of the overlap region, as indicated by the arrow. 




R c (A/m) 

FIG. 3. The total free energy of the binary mixture, F, 
when the large (0.237-/im) sphere was at the vesicle wall. F 
is plotted in units of feT as a function of the wall curvature 
radius, Re- The symbols represent measurements and the 
line represents the results of the calculation described in the 
text. (Inset) F/ksT vs. position, s, along the perimeter of 
the vesicle. The origin (s=0) is located at the lower left of 
the vesicle shown in Fig. 1(c). 




FIG. 4. Cartoon demonstrating how hard spheres can in- 
duce a change in shape of a vesicle, (a) When the membrane 
is flat, the excluded- volume overlap is minimized and the cur- 
vature energy is zero, (b) When the vesicle envelops the large 
sphere, the curvature energy increases, but so does the ex- 
cluded- volume overlap (hatched region). 
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